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' Abstract 

o ■ 

' I review some recent developments in the QCD description of elas- 

tic vector meson production, focusing on issues like the meson wave 
function, the gluon density in the proton, factorisation, the use of 
parton-hadron duality and the nonperturbative vacuum. 



Oh! 

: 

■ 1 Introduction 

5^ ! I will discuss some recent developments in the description of elastic vector 

meson production at high energy within QCD, restricting myself to the cases 
where either the photon virtuality or the vector meson mass My is large 
while the squared momentum transfer t from the proton is small. I will high- 
light progress and problems in the theory of these processes, concerning issues 
such as the meson wave function, the gluon density in the proton and factori- 
sation. I shall also spend some time on new approaches, namely the use of 
parton-hadron duality and a nonperturbative description of the scatter- 
ing 0. My aim is not to attempt a detailed comparison of theory predictions 
with data but rather to indicate where and why the predictions of various 
authors are different. I will not have time to speak about phenomenological 
models based on Regge theory. 

*Talk given at the Ringberg Workshop on New Trends in HERA Physics, Ringberg 
Castle, Germany, 25-30 May 1997 

**Unite propre 14 du Centre National de la Recherche Scientifique. 



1 




Figure 1: One of the Feynman diagrams for ^y* + p ^ V + p with a vector 
meson V = p, p', u, 0, J/ip, etc. More diagrams are obtained by attaching 
the gluons to the quark loop in different ways. 

2 Simple approach 

The common picture of most QCD models for elastic vector meson production 
at high energy is shown in Fig. |l|, where I also define the kinematics. At 
large 7*p cm. energy W the dominating exchange is the pomeron, which is 
described by the exchange of two gluons coupling on one side to the proton 
and on the other to a quark loop to which are attached the virtual photon 
and the vector meson. The vector meson couples to the quark line of the 
loop by its wave function il^lz, kx)- Here z and fc-r are defined by a Sudakov 
decomposition for the quark momentum in the meson, k = zq' + Cp + ^t- 

As a baseline let me present a very simple version of this model, recently 
used by Cudell and Royen 0]. For the meson it takes a nonrelativistic, 
constituent quark description, namely a wave function peaked at z = ^ 
and kT = 0, so that the quark and antiquark with mass My/2 share the 
four-momentum of the meson equally. It makes the approximation that the 
exchanged gluons do not interact and couple directly to the three constituent 
quarks of the proton, following the work of Low, Nussinov, Gunion and 
Soper §. 

With these assumptions Cudell and Royen find the Q^- and t-dependence 
of the cross sections for p, and J/ip production in fair agreement with 
data from HERA, EMC and NMC, except for problems with the fixed tar- 
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get J/ip data. Let me remark that the cross sections for longitudinal and 
transverse photons respectively behave like ai ~ l/Q^ ~ fo^' 

^ My, while for nonasymptotic the behaviour in this variable is more 
complicated, even in this simple model where a fixed strong coupling was 
taken. One finds however a ratio oi^joT = Q"^ /My for all > 0, which is 
in disagreement with the data — a point I will come back to in Sec. 4. As 
the approximation of noninteracting gluons does not allow to describe the 
energy dependence of the cross section the authors introduce an overall phe- 
nomenological factor in ai and cjt, and find consistency with the data if this 
factor is allowed to depend on W, but not on or My. Little change in the 
results is obtained when the propagators of the exchanged gluons are taken 
as nonperturbative, as was done in earlier work of Donnachie, Landshoff and 
Cudell Finally, the authors observe that their model agrees with the 
data down to = even for the light mesons p and (p. Notice that in this 
case one has no large virtuality in the quark loop to justify the perturbative 
treatment of the quarks and their coupling to the photon. 

3 Sophistication of the model 

Having seen that even the simplest form of the model depicted in Fig. ^ does 
not fare badly in describing data we will now take a look at several of its 
refinements. 

3.1 Meson wave function 

The meson wave function i/j{z, kr) turns out to be a major source of un- 
certainty in the predictions of the model. It determines which virtualities 
dominate the integration over the quark loop and in particular influences the 
overall normalisation of the cross section, its Q^-dependence and the pro- 
duction ratios for different mesons. Various choices have been made in the 
literature ^ |^, ^, |^, [1^, [TT| , sometimes leading to appreciable differences 
in numerical results. For J/ip photoproduction, as an example, Frankfurt 
et al. [0] find a suppression of the cross section due to finite fcr in the wave 
function which is significantly stronger than the one estimated by Ryskin et 
al. 0. For p electroproduction a wide range of suppression factors has been 
obtained by Frankfurt et al. according to the choice of ilj{z,kT), with a 
particular sensitivity to its large-fcy tail. 

On the theory side it is known that as the renormalisation scale for the 
wave function tends to infinity, corresponding to very large in the pro- 
cess, the /cT-integrated wave function / dh'^iplz, kx) behaves like z{l — z) 
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due to QCD evolution, both for longitudinal and transverse meson polari- 
sation. This is very different from the 5{z — |) in the nonrelativistic wave 
function, but the difficult question is of course to assess how far one is from 
the asymptotic regime for the values of one has in experiment. Inves- 
tigations of the wave function using the operator product expansion and 
QCD sum rules have been made by Ball and Braun [jl2| and Halperin and 
Zhitnitski UTTl. The latter find that at z — ^ and z ^ 1 the wave func- 



tion should only depend on ^(Yz^, which in particular excludes a factorising 
ansatz ipi^z^ kx) = 4>{z) -xikr)- They also argue that from a conceptual point 
of view there should be no perturbative tail like l/k"^ in the wave function 
and that the corresponding corrections to the meson-quark vertex should be 
explicitly treated as a^-corrections to the leading order result. 

In the operator product expansion framework the effects of finite kj- in 
il){z,kT) are of higher twist. It is important to realise that there are other 
higher twist contributions from diagrams where one or more gluons enter 
the wave function blob in Fig. 0. They correspond to higher fock states of 
the meson such as qqg and qqgg^ which have respectively been considered 
by Halperin and Hoodbhoy |TD|. Moreover, taking only into account a 



qq wave function ip{z,kT) is not gauge invariant. To see this consider that 
a finite fcr involves the transverse component of the operator between 
quark fields in a matrix element. To recover gauge invariance one must pass 
to the covariant derivative, + A^, whose transverse component involves 
physical gluon degrees of freedom. 

Finally let me mention a problem raised by Frankfurt et al. concern- 
ing the J/ip wave function. Calculating the transition amplitude of a. J/ip 
to a timelike photon via a charm quark loop with wave functions 'ip{z^kT) 
obtained in various nonrelativistic potential models they find that the region 
of kf larger than the charm mass contributes up to 30% of the loop integral, 
thus calling into question the consistency of a nonrelativistic description. 
Note that the transition J/?/' — > 7* ^ e^e~ is commonly used to fix the 
normalisation of the wave function from the semileptonic J/ip decay width. 



3.2 Quark mass 

Which quark mass should be chosen in the loop of Fig. |] is another source of 
uncertainty. Surprisingly this even holds for heavy mesons. Ryskin et al. |^ 
point out that the cross section scales like the eighth power of Mj/^/{2mc), 
resulting in huge changes of its overall normalisation even if one allows for 
a modest change in the charm mass m^. Incidentally, Frankfurt et al. 
advocate to use the running quark mass rather than the pole or constituent 
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mass. 

In the case of hght quarks both the choice of current and constituent 
masses has been made in the literature |, 0. Considering the limit 

of zero quark mass rriq reveals that light meson production from transverse 
photons is a higher twist effect the cross section ax is proportional to 
in the coUinear approximation for the quarks, i.e. if one sets their transverse 
momentum kx to zero in the calculation of the loop and uses the /cT-integrated 
meson wave function, corresponding to what one does when using collinear 
parton densities in DIS. At least one of the quantities or kx must be 
nonzero to obtain a nonvanishing <Tx from diagrams as shown in Fig. |^. 



3.3 Off-diagonal gluon density 

We now focus on the lower blob in Fig. |I], the coupling of the two gluons 
to the proton, including their interaction. If the outgoing proton had the 
same momentum as the incoming one this blob would be given by the gluon 
density in the proton. In our reaction this is however not the case: even if 
the proton has zero transverse momentum after the scattering it has lost a 
fraction x = {My + Q'^)/{W'^ + Q"^) of its longitudinal momentum, which is 
necessary to make a timelike vector meson out of a spacelike photon. Thus 
the longitudinal momentum fractions Xi and X2 of the gluons with respect 
to p (cf. Fig. ID) are not equal, their difference being xi — X2 = x. One can 
argue that in the leading InQ^ ■ln(l/a;) approximation it is indeed the gluon 



density g{x) that describes the blob |jT^, |T^, as the typical values of Xi and 
X2 in the loop integration are of order x and to leading ln(l/a;) all such values 
are equivalent. 

The most prominent phenomenological consequence of this description is 
the strong rise of the cross section in W when either of the scales or My 
is hard. At this point I should remark that only the imaginary part of the 
7*p — > Vp amplitude is actually calculated from the Feynman diagrams of 
Fig. |l|, cutting them in the s-channel. The contribution of the real part to 
the cross section is negligible if the W^-dependence is weak as in the case of 
soft pomeron exchange, otherwise it can be accounted for in an approximate 
way [H i. 

Beyond leading ln(l/x) the quantity describing the blob in Fig. |1] is cer- 
tainly different from the gluon distribution that contributes to inclusive DIS, 
and there has been much theoretical interest recently in such so-called "off- 
diagonal parton densities" . Strictly speaking they are not "densities" by the 
way, since they correspond to matrix elements of parton fields between dif- 
ferent proton states and do not have a probability interpretation. Let me 
note that their evolution equations differ from the DGLAP equations by their 
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splitting functions |T6|, |I7| . Some studies have been performed [|T^, [T^ to es- 
timate the difference between the off-diagonal and the usual, diagonal gluon 
densities, finding no big difference at small x, but much remains to be done 
in this field. An obvious remark is that the dependence of the off-diagonal 
gluon density on t is specific of its asymmetric kinematics and cannot be 
predicted by approximating it with the diagonal density measured in other 
processes. 

In the same way as in DIS at small x one can also go beyond the leading 
InQ^ approximation and instead of g{x) consider the unintegrated gluon 
density f{x, It) with a finite transverse momentum It of the gluons 0. 

The question of which factorisation scale is appropriate in the gluon den- 
sity is numerically very important since we know from DIS that parton dis- 
tributions at small x change rapidly with this scale. Generically it is set 
by the typical virtualities in the quark loop, i.e. by and the quark mass, 
but various concrete choices have been advocated in the literature: for J/tp 
production the scale ^(Q^ + Mj^^) used by Ryskin et al. ^ is smaller than 
that of Frankfurt et al. , whereas Nemchik et al. choose + My times 
a factor between 0.07 and 0.2 depending on the meson and its polarisation. 
Let me however emphasise that any choice of factorisation scale can only be 
an educated guess of which value will make higher order corrections small. 
As in other reactions such as for instance jet production experiment should 
feel free to try different scales const - {Q^ + My), and also const - for light 
mesons. 



3.4 Factorisation theorem 

For the asymptotic limit where is much larger than all masses in the 



process a factorisation theorem has been proven by Collins et al. |T8| within 
perturbative QCD. By an analysis of Feynman diagrams and power counting 
arguments they find that to all orders in the amplitude for 7*p — >■ Vp 
factorises into a hard scattering part, a collinear, off-diagonal quark or gluon 
distribution in the proton, and a collinear qq wave function for the vector 
meson (see Fig. |]). The corresponding amplitude has a power behaviour like 
1/Q. 

This power behaviour and factorisation property holds only for a longitu- 
dinal photon. If the 7* is transverse the authors find that the amplitude only 
behaves like and that there is no factorisation like in Fig. |^: soft, "wee" 

partons can bypass the hard scattering and directly go from the proton blob 
A to the meson blob B. The power behaviour in Q confirms from a differ- 
ent point of view the statement that light meson production from transverse 
photons is a higher twist process. 
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Figure 2: Factorisation of the amplitude for 7*p — > Vp with a longitudinal 
photon into a hard scattering part H, an off-diagonal quark or gluon density 
A and a meson wave function B. 

It is worth noting that these findings are not restricted to the small-x 
region but valid for general x, where the quark distribution in the proton is 
relevant in addition to that of the gluons. 

Finally it is clear that the direct phenomenological application of these 
results requires "sufficiently large" values of Q^, and we have already seen 
that higher twist effects such as the transverse momentum in the meson 
wave function can be numerically important in the HERA regime. Also one 
should keep in mind that the Q^-behaviour obtained by power counting can 
be strongly masked by logarithmic scaling violations which are strong at 
small X. 



3.5 Rescattering effects 



Let me just briefly mention that several papers p, ||, |1^ have investigated 



the effects of rescattering of the gg-pair on the proton, in other words shad- 
owing corrections or multiple-pomeron exchange, and found that they can be 
important in kinematic situations accessible at HERA. An observable that 
is naturally sensitive to such effects is the t-dependence of the cross section, 
which has been investigated in detail by Gotsman et al. [pJ[| . 

In this context one may remark that an exponential parametrisation 
da/dt oc exp{bt) at small and moderate values of t is no more than a fit to 
a simple function: from theory one does not expect such a behaviour to be 
exact. Remember for instance that ordinary elastic form factors of hadrons 
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are not exponentials. One should be aware that just comparing an experi- 
mentally fitted slope parameter b with the logarithmic slope of da / dt a.t t = 
calculated in theory can be misleading. To compare the full t-dependence is 
of course the most thorough way to proceed, but if one looks for convenient 
"handy" parameters other choices than b might be useful, such as the mean 
value {t) proposed by Cudell and Royen 0. 

4 The ratio R = (Jl/ (Jt 

As I have already emphasised the physics of light meson production is quite 
different with transverse and longitudinal photon polarisation. It turns out 
that (Jr is more sensitive to small quark virtualities and thus to infrared 
physics than its counterpart for longitudinal photons H, ||. This can 
be understood from the properties of the transition 7* qq: while for 
a longitudinal photon configurations are preferred where the longitudinal 
momenta of quark and antiquark are comparable, a transverse photon likes to 
split into a gg-pair where the quark or antiquark carries only a small fraction 
of the photon momentum and is soft. This is the aligned jet configuration, 
which is of prime importance for the physics of diffraction. The interaction 
with the gluons hardly changes the longitudinal momenta of q and g, so 
that this configuration corresponds to small 2; or 1 — 2; in the meson wave 
function. To which extent (Tt is dominated by the soft region thus depends 
on how strongly small z or 1 — z are suppressed in ip{z, k^). 

In the simple model presented in Sec. 2, where a meson wave function 
peaked at 2; = | is used, the ratio R comes out as Q"^ /My, which is far too 
big to describe the data for in the range of some GeV^. In the model 
of Nemchik et al. |p the greater infrared sensitivity of ctt is reflected in a 
lower factorisation scale of the gluon density than in ctl, and the authors find 
that R grows less fast than linearly in and that its rate of growth in 
depends on x. 

4.1 Using parton-hadron duality 

To circumvent the uncertainties associated with the meson wave function in 
this delicate context, Martin et al. invoke parton-hadron duality: to obtain 
the cross section for p-production they calculate the diffractive production 
of an open gg-pair by two-gluon exchange, project out the appropriate gg 
partial wave to make a vector meson of given helicity, and integrate the 
corresponding cross section over the gg invariant mass Mqq in an interval 
around Mp. The authors stress that in this mass range partonic final states 



8 



other than qq are heavily suppressed and that the dominating hadronic final 
state is a pair of pions. In the ratio R several uncertainties are expected to 
cancel to a large extent, e.g. uncertainties about the appropriate interval of 
Mqq and about radiative corrections (the authors estimate that there should 
be a large i^'-factor). As in the calculation by Nemchik et al. ay is 
dominated by lower quark virtualities than cxi, now in the context of open qq- 
production, and through the different factorisation scales in the gluon density 
the authors obtain a ratio R that grows less than linearly in and is in fair 
agreement with HERA data for R in the range 5 GeV^ < g2 < 20 GeVl 
Again this ratio is predicted to depend on x. 



5 A nonperturbative model 

To finish let me present an approach from a perspective of nonperturbative 
physics. The model of Dosch et al. starts from the high-energy scattering 
of a quark or antiquark, which is described as the scattering in an external 
gluon field pO], like in the semiclassical model of diffraction of Buchmiiller 



and Hebecker |21[]. This gluon field is then averaged over in the sense of a 
path integral, according to the stochastic model of the nonperturbative QCD 
vacuum High-energy scattering thus becomes related with nonpertur- 
bative parameters such as the SVZ gluon condensate. 

The nonabelian character of the gluon field is crucial in this model: not 
single quarks or antiquarks are scattered but rather colour singlet configura- 
tions of qq or qqq, their transverse separation having a strong infiuence on the 
strength of the scattering. The physical picture for vector meson production 
is then as shown in Fig. ^ For the three-quark wave function of the proton 
and the qq wave function of the vector meson an ansatz has to be made, 
whereas the splitting of the virtual photon into qq is calculated perturba- 
tively. As in the simple two-gluon exchange model discussed in Sec. 2 the 
IV-dependence of the cross section cannot be predicted in the present form 
of the model. 

Its authors find fair agreement of their results with EMC and NMC data 
for Q"^ between 2 and 10 GeV^ and W between 10 and 20 GeV, looking at 
observables such as the Q^- and t-dependence of the cross section and its 
normalisation, except that their cross section for the is about a factor 2 
above the data. They also describe the ratio R for p-production measured 
by NMC and the t-dependence of J /if) photoproduction. They conclude by 
emphasising the need to incorporate nonperturbative aspects of the photon 
wave function for small and My, and perturbative gluon contributions 
which they expect to become important as W increases from fixed target to 
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Figure 3: The process 7*p Vp in the model of Dosch et al. |@]. Quarks 
and antiquarks are scattered in the gluon field of the nonperturbative QCD 
vacuum. 

HERA energies. 

6 Conclusions 

The model depicted in Fig. [l| is a good candidate for the description of elastic 
vector meson production at large or meson mass in QCD. Even a very 
simple version of this model as presented in Sec. 2 can describe quantitative 
features of the data. 

Progress has been made in the theory of meson wave functions and of 
the off-diagonal gluon distribution in the proton. For the asymptotic regime 
of very large a factorisation theorem has been worked out which further 
elucidates the different physics of ai and ar- 

On the phenomenological side various difficulties or uncertainties present 
themselves if one is to make quantitative predictions in the nonasymptotic 
regime, in particular for light meson production from transverse photons, 
or if in turn one aims to extract information from the data on the meson 
wave function or the gluon distribution. An important task is to identify 
observables which are sensitive to only a few theoretical effects. 

Finally, as I showed in Sec. 4.1 and 5, there are promising alternative 
approaches, which highlight the connection of this reaction with open qq- 
production, and with nonperturbative physics and the QCD vacuum. 
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